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We show that the critical behaviour of two- and three-dimensional frustrated magnets cannot 
reliably be described from the known five- and six-loops perturbative renormalization group results. 
Our conclusions are based on a careful re-analysis of the resummed perturbative series obtained 
within the zero momentum massive scheme. In three dimensions, the critical exponents for XY 
and Heisenberg spins display strong dependences on the parameters of the resummation procedure 
f) |' and on the loop order. This behaviour strongly suggests that the fixed points found are in fact 

spurious. In two dimensions, we find, as in the O(N) case, that there is apparent convergence of 
the critical exponents but towards erroneous values. As a consequence, the interesting question 
of the description of the crossover/transition induced by Z2 topological defects in two-dimensional 
frustrated Heisenberg spins remains open. 
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d • I. INTRODUCTION 
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After more than thirty years of intensive studies, the critical behaviour of frustrated magnets is still 
controversial (see [l| and references therein) . At the root of the problem is the competition between 
the interactions among neighboring spins that gives rise to a canted ground state, and thus, to a 
1-^ , symmetry breaking scheme where the rotational group is fully broken. This is for instance the case 
r"j ' in the paradigmatic example of frustrated magnets, the Stacked Triangular Antiferromagnets (STA), 
O ' where the three spins on an elementary cell display a planar 120° structure in the ground state. As a 
consequence the order parameter is a matrix instead of a vector (a 5*0(3) matrix for Heisenberg spins 
and a 2 x 2 matrix for XY spins) and the critical properties are therefore entirely different from those 
— , of unfrustrated systems. 

For instance in dimension d = 2, the first homotopy group of SO(3) being non trivial - Ui(SO(3)) = 
Z2 - one expects for Heisenberg spins a deconfinement of topological excitations that could give 
rise to a Kosterlitz-Thouless(KT)-like transition @,J3 or, at least, to a crossover behaviour - see 
below. Numerous experimental [5|-|9( and numerical [2j, Il0l4l7l | studies have indeed shown indications 
of a nontrivial phenomenon occuring at finite temperature. For XY spins the order parameter space 
is given by SO(2) x Z2. In this case coexist Ising degrees of freedom, topological excitations - 
Hi(SO(2)) = Z - and spin-waves. A very debated issue has been the nature of the phase transition(s) 
occuring in this system as the temperature is varied: either two separate Ising and KT transitions or 
a unique one (see [l8| and references therein). 

In d = 3 the question of the criticality has been extremely controversial (see [l|). On the one hand, 
many experiments display scaling behaviours for XY and Heisenberg spins with critical exponents 
differing from those of the O(N) universality class (see [l[ for a review). On the other hand, other 
experiments as well as extensive Monte Carlo simulations performed on STA or on similar models have 
exhibited weak first order behaviour [l9l - l27j . Two main explanations have been proposed to describe 
these contradicting results. 

The first one is based on a perturbative renormalization group (RG) approach performed at fixed 
dimension (FD) either within the minimal-substraction (MS) scheme without e-expansion [28| or 
within the zero momentum massive scheme [29|-|32j, at five- and six- loop order respectively. Within 
these approaches, stable RG fixed points were found for N = 2 and N = 3 leading to the prediction of 
second order phase transitions in d — 3 for frustrated magnets. Note that, within these FD approaches, 
one also finds a fixed point in d = 2 with nontrivial critical exponents in the N = 2 and N — 3 cases 
[3ll I33L HiJ . This fact has led to the hypothesis of a Kosterlitz-Thouless-like behaviour induced by Z2 
topological defects for Heisenberg spins [331 ] . 

The second explanation is based on both the e — 4 — d (or pseudo-e)-expansion [35U37| and the 
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nonperturbative renormalization group (NPRG) approaches [H. l38M4ll | . In these approaches, one finds 
that there exists, within the (d, N) plane, a line N c (d) that separates a second order region for N > N c 
from a first order region for N < N c . Within the e-expansion, one finds N c (d = 3) ~ 5.3 J37j. pseudo-e 

.1 [Jsc 



expansion gives N c (d — 3) ~ 6.23 36] and within the NPRG approach N c (d — 3) ~ 5.1 1] so that the 
transition for N = 2 and N = 3 are predicted to be of the first order. A thorough analysis performed 
within the NPRG approach [l|, [H, Hoi ] has shown that even if there is, strictly speaking, no fixed 
point below N c (d = 3) the RG flow is very slow for N = 2 and 3 in a whole region of the coupling 
constant space so that there is pseudo-scaling without universality on a large range of temperature, 
in agreement with the numerical and experimental data. 

Although the NPRG approach very likely explains the whole body of known data (see [l[ for details), 
several points of the physics of frustrated magnets remain controversial. The first one is that, although 
the occurence of (weak) first order transitions is by now well established in several three-dimensional 
systems this does not imply that all systems, sharing the same order parameter and the 

same symmetries, undergo first order phase transitions (see for instance [H[ for a recent numerical 
computation where a second order phase transition is observed). In other words, the very existence 
of a parameter domain where the frustrated systems would undergo second order transitions is still 
debated. A second point that should be understood is the origin of the discrepancy between the 
two scenarii above and, in particular, between the results obtained within the different perturbative 
schemes: e (or pseudo-e)-expansion on the one hand and the FD approaches on the other hand. A 
last and important question is the nature of the transition that occurs in d = 2 for Heisenberg spins: 
phase transition or simple cross-over behaviour between a low-temperature - spin-wave - phase and 
a high-temperature phase with both spin-waves and vortices ? 

It is clear that answering to these questions amounts to answering to the question of the existence of 
a genuine attractive fixed point in the RG flow of frustrated magnets in d = 3 and 2 for N — 2 and 3. 
In the perturbative framework, this is essentially equivalent to proving (or disproving) the reliability 
and convergence properties of the resummation procedures necessary to obtain sensible results out 
of the perturbative series. This work has been initiated in d = 3 in a previous publication where 
the five-loop perturbative series obtained in the MS scheme without e-expansion have been carefully 
reexamined |42j ] . Studying (i) the convergence properties of the critical exponents with the order L 
of the expansion (number of loops) and with respect to the variations of the parameters involved 
in the resummation procedure, (ii) the properties of the fixed point coordinates (u*(d,N),v*(d,N)) 
considered as functions of d and N and (iii) the continuation of the fixed point found in d — 3 for 
N = 2 and N — 3 up to d = 4, the authors of 42j have provided strong arguments in favor of the 
spurious character of the fixed points found in d = 3, i.e. that they are artefacts of the perturbative 
expansion in the MS scheme. 

In the present paper, we extend the previous analysis to the series obtained in the zero momentum 
massive scheme in d = 3 at six loops (30| and in d = 2 at five loops [lij]. In d = 3, we apply the 
criteria used in [42|, |43[ — Principle of Minimal Sensitivity (PMS) and Principle of Fastest Apparent 
Convergence (PFAC) — and confirm that the fixed points found in d = 3 for N = 2 and N = 3 are 
most likely spurious. In d = 2, the situation is more delicate. We recall that, already for the (non- 
frustrated) O(N) models, the critical exponents found from the (<j> 2 ) 2 field theory are quantitati vely 
wrong although apparently converged for all N. This striking phenomenon, already mentioned in [44|, 
relies on the presence of non-analytic contributions to the /3-function at the fixed point [45| . The same 
kind of problem has been mentioned in the case of frustrated magnets [U H3| but it was assumed to 
leave unaffected the qualitative predictions, in particular, the existence of a non-trivial fixed point in 
the Heisenberg case. We show here, on the contrary, that the phenomenon of apparent convergence 
towards erroneous values, together with the presence of instabilities of some critical exponents with 
respect to the resummation parameters, leads to seriously question the conclusions drawn in the past 
as for the critical behaviour of these systems. 

Our study altogether shows clear evidences that the FD perturbative approaches to three- 
dimensional frustrated magnets are not reliable, at least at the orders studied, and that there is 
no convincing evidence of a genuine phase transition induced by vortices in two-dimensional Heisen- 
berg spin systems. Finally, for the same reasons as in the Heisenberg case, we show that the behaviour 
of XY spins in d = 2 cannot be elucidated from the five-loop perturbative data. 

The paper is organized as follows. In Section II, we study in detail the O(N) case in d = 2 and 
d = 3. In d = 3, this allows us to illustrate on a well-known example the kind of stability (resp. 
instability) properties expected for a genuine (resp. spurious) fixed point. In d = 2, this allows us to 
illustrate the fact that there can be fast apparent convergence of the critical exponents but towards 
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erroneous values due to nonanalytic contributions. In Section III, frustrated magnets are studied in 
d = 3 and in d = 2. In d = 3, we confirm the spurious character of the fixed points found for N = 2 
and N — 3. In d — 2, we show the unreliability of the conclusions - phase transition controlled by a 
fixed point - deduced from the results obtained perturbatively at five loops in the N = 3 case. We 
then analyze the N = 2 case and reach the same conclusions as in the N — 3 case. 



II. THE O(N) MODELS IN TWO AND THREE DIMENSIONS 

In the following, we study the convergence of the resummed perturbative series obtained for the 
frustrated models. Since we need to determine criteria to decide whether the perturbative results are 
(or are not) converged, we illustrate briefly how convergence of the resummed series shows up for 
the O(N) models in d — 3. We show that the behavior of the correction to scaling exponent w as a 
function either of the loop order or of the resummation parameters is a good indicator of the numerical 
convergence of the perturbative results. The exponent ry, when available, is also studied. By analyzing 
the two-dimensional O(N) case we also show that, contrarily to common belief, the five-loop results 
for the critical exponents are not converged. The reason of this behaviour is however rather subtle 
since there is, in fact, apparent convergence but towards erroneous values, a phenomenon that we call 
anomalous apparent convergence. 



A. Resummation procedure 

As well known, the perturbative series obtained in the O(N) models for the /3 function describing 
the running of the coupling constant with the scale are not convergent (46l . |47| . They are asymptotic 
series which, in the case of the zero momentum massive scheme, are Borel summable [48]. Powerful 
resummation methods have been used in the past that, thanks to the knowledge of the behavior of 
the series at large order, lead to converged and accurate results (see [U S^] for reviews). We recall 
in the following the kind of resummation procedure that we use throughout this article. 

Let us consider a series 



f(u)=J2a n u n (1) 

n 

where the coefficients a n are supposed to grow as n\. 
The Borel-Leroy sum associated with f(u) is given by: 

B( u ) = Y u n (2) 

n 1 J 

where 6 is a parameter whose meaning will become clear later. 

The resulting series is now supposed to converge, in the complex plane, inside a circle of radius 1/a, 
where u = — 1/a is the singularity of B(u) closest to the origin. Then, using this definition as well as 
r[n + 6+1] = J °° t n+b e^ t dt, one can rewrite 

/>oo 

/(„) = y —^1 u n / dt e-* t n+b . (3) 



T\n + b- 

Interchanging summation and integration, one can now define the Borel transform of / as: 

poo 

f B (u) = dt e -t t b B(ut) . (4) 



In order to perform the integral Q on the whole real positive semi-axis one has to find an analytic 
continuation of B(u). Several methods can be used, Pade approximants constitute one possibility 
[50M52I ] . However, it is generally believed that the use of a conformal mapping [53|, [HJ is more efficient 
since it makes use of the convergence properties of the Borel sum. Under the assumption that all 
the singularities of B{u) lie on the negative real axis and that the Borel-Leroy sum is analytic in the 
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whole complex plane except for the cut extending from — 1/a to — oo, one can perform the change of 
variable: 



uj(u) = 



\/\ + au — 1 



VI + au+ 1 



u(us) — 



a(l-w) 2 



(5) 



that maps the complex w-plane cut from u = —1/a to — oo onto the unit circle in the w-plane such 
that the singularities of B(u) lying on the negative axis now lie on the boundary of the circle \w\ = 1. 
The resulting expression B(u(lo)) has a convergent Taylor expansion within the unit circle \u)\ < 1 
and can be rewritten: 



fl(u(w)) = 5^d n (a,&) Mu)} 7 



(6) 



where the coefficients d n (a, b) are computed so that the re-expansion of the r.h.s. of ^ in powers of 
u coincides with that of ([1]). One obtains through © an analytic continuation of B(u) in the whole 
u cut-plane so that a resummed expression of the series / can be written: 

/>oo 

/k(«) = Y] Ma, b) / dt e-* t b [u>(ut)] n . (7) 
In practice it is interesting to generalize the expression (J?} by introducing [55| the expression 



fR( u ) = y"^ n (q,a, b) 



dt e - * t b 



Hut)] n 

[1-L0(ut)] a 



(8) 



whose meaning will be explained just below. 

If an infinite number of terms of the series //j(u) were known, expression would be independent 
of the parameters a, b and a. However when only a finite number of terms are known, fn(u) acquires 
a dependence on them. In principle, the parameters a and b are fixed by the large order behavior of 
the series: 



(-a lo ) n n\n b 



(9) 



which leads to a = a\ and b > b\ a + 3/2 [54| where a\ a and &i denote the large-order value of a and 
6. As for a, it is determined by the strong coupling behavior of the initial series: 



(10) 



which can be imposed at any order of the expansion by choosing a = ao. The common assumption is 
that the above choice of a, b and a improves the convergence of the resummation procedure since it 
encodes exact results. 

Let us however emphasize that, often, only a is known and that the other parameters, a and 6, 
must be considered either as free (as for instance in (28|) or variational (as for instance in |42l. |43| 
where a is determined by optimizing the apparent convergence of the scries). In any case, the choice 
of value of a, a and b must be validated a posteriori. 



B. O(N) models in three dimensions and principles of convergence 



The dependence of the critical exponents upon the parameters a, b and a is an indicator of the 
(non-) convergence of the perturbative series. Indeed, in principle, any converged physical quantity Q 
should be independent of these parameters. However, in practice, at a given order L of approximation 
(loop order), all physical quantities depend (artificially) on them: Q —> Q^(a, 6, a). Even if a is fixed 
at the value obtained from the large order behavior, all physical quantities remain dependent upon b 
and a at finite order. We consider that the optimal result for Q at order L corresponds to the values 
(bjp t , a^l) of (o, a) for which Q depends most weakly on b and a, i.e. for which it is stationary: 



Qopt ~ H^ipti a opt) 



with 



dQ<- L \b, a) 



db 



dQl L >{b,a) 
da 







(11) 
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where, of course, and a^l are functions of the order L. 

The validity of this procedure, known as the "Principle of Minimal Sensitivity" (PMS), requires 

that there is a unique pair (&opt> a opt) sucn tnat lii stationary. This is generically not the case: 
several stationary points are often found. A second principle allows us to "optimize" the results 
even in the case where there are several "optimal" values of b and a at a given order L: this is the 
so-called "Principle of Fastest Apparent Convergence" (PFAC). The idea underlying this principle is 
that when the numerical value of is almost converged (that is L is sufficiently large to achieve a 
prescribed accuracy) then the next order of approximation must consist only in small change of this 
value: Q( L+1 ) ~ Q^ L \ Thus, the preferred values of b and a should be the ones for which the difference 
between two successive orders Q^ L+1 ^ [b^ L+1 \ o/ L+1 ') — Q^(b^ L \ a^ L >) is minimal. In practice, the 
two principles should be used together for consistency and, if there are several solutions to Eq. fTTj) at 

order L and/or L + 1, one should choose the couples (&opt' a opt) an d (&opt ! a opt ) f° r which the 

stationary values Q^(&opt> a opt) an d ( 9^ i+1 H^opt" 1 \ v-opt^) are tne closest, that is for which there is 
fastest apparent convergence. These principles have been developed and used in [42|, [43|, [54[ , see also 

M- 

Nice examples where these two principles work very well and indeed lead to optimized values of the 
critical exponents are the O(N) models in d — 3 computed perturbatively at four-, five- and six- loop 
orders (within the zero momentum massive scheme). The series for the /3- function of the coupling 
constant are resummed thanks to a conformal Borel transform. Subsequently, one obtains the fixed 
point coordinate u*, its stability being defined by the correction to scaling exponent ui: 



d(i{u) 



On 



(12) 



A positive value of ui (or a positive real part if it is complex) corresponds to a stable fixed point. 
We show in FigQ]the exponent ui of the 0(4) model in d = 3 as a function of the parameter b for 
the values of a for which stationarity is found for both b and a. As expected, the dependence of uj 
upon the resummation parameters becomes smaller as the order in the loop expansion increases as 
illustrated by the curves ui(b) that flatten between four and six loops, see FigfT] At this order one 
finds uj ~ 0.783. As an indicator of the quality of the convergence we give the difference between the 
fifth and the sixth order for the exponent ui: ui(L — 6) — ui(L — 5) ~ 2 . 10 -4 . Note that this case 
also illustrates the situation where - at six-loop order - several stationary points occur and where the 
PFAC allows us to select a single solution, see FigfT] Our results are comparable with six-loop results 
obtained by Guida and Zinn- Justin [13] for N = 4: uj = 0.774 ± 0.020 (in d = 3), uj = 0.795 ± 0.030 
(within the e-expansion) . 




FIG. 1: The exponent cj of the three-dimensional 0(4) model as a function of the resummation parameter b 
at four-, five- and six-loop orders. The dot on each curve corresponds to a stationary value of uj = ui{a,b) in 
both a and b directions with a fixed to its large-order value a\ — 0.1108. One has: (a op t, b opt ) = (3.2, 8.5) at 
four loops, (a op t,6o P t) = (5.2,9.5) at five loops and (a pt,&o P t) = (5, 13) at six loops. 

It is remarkable that the same study performed on other critical exponents, other values of N and 
even with other perturbative series (obtained from the MS scheme for instance) always leads to the 
same kind of results with values of the exponents found that are very close to the best known values 
obtained from Monte Carlo simulations. This proves that the above methodology is indeed efficient. 

Finally note that the same argument can also be applied to the determination of an optimal value 
of a, a opt , from the PMS applied to this parameter: if there is convergence of the resummed series, 
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we expect that a op t almost coincides with the value determined by the large order analysis, Eq.©, 
a opt — a\ . The difference between these quantities is a measure of the convergence level of the series. 
We show in Figf2]on the example of the 0(4) model in d = 3 that, as expected, the value a op t is very 
close to a\ and the difference between w(a pt) and u>(a\ ) is extremely small. 

0.784 
0.782 
0.780 
0.778 

0.10 0.11 0.12 0.13 




FIG. 2: The exponent ui of the three-dimensional 0(4) model as a function of the resummation parameter a 
at four, five and six loops. The vertical line corresponds to a — a\ — 0.1108. The values chosen for a and b 
are such that uj is stationary w.r.t. a and b when a = a\ (see FigQ}. Note that w(a op t) — oj(a\ ). 



The criteria given above are of crucial importance, especially when considering FD approaches. 
Indeed, generically the (non-resummed) series obtained at L loops for the /3-function are polynomials 
of order L + 1 in the coupling constant u. Thus, the fixed point equation /3(u*) = admits L + 1 
roots u* that are either real or complex. Contrary to what occurs within the e-expansion, where the 
coupling constant is by definition of order e, in the FD approach no real root can be a priori selected 
or, reciprocally, discarded. As a result, the generic situation is that the number of fixed points as 
well as their stability vary with the order L: at a given order, there can exist several real and stable 
fixed points or none instead of a single one. In principle, the resummation procedure allows both to 
restore the non-trivial Wilson-Fischer fixed point and to eradicate the non-physical, spurious, roots. 
In particular, we expect that spurious solutions should satisfy neither the PMS nor the PFAC criteria. 
We show in Fig|3]on the example of the 3d 0(A) model that there exists, beside the usual Wilson- 
Fisher stable fixed point, a spurious (unstable) fixed point. As expected, the exponents computed 
at this fixed point are very unstable w.r.t. variations of the resummation parameters a (the same 
behaviour occurs when variations of b are considered), a behaviour which, according to our criteria, 



is sufficient to discard it [77 1. 




FIG. 3: The exponents uj and 77 as functions of a at the spurious fixed point of the 0(4) model in d — 3 at 
four and six loops. The vertical line corresponds to a — a\ — 0.1108. The values considered for a and b are 
respectively equal to 6 and 4. Other values gives similar results. The exponent to is negative since the spurious 
fixed point is repulsive. 



A potential difficulty with this procedure is that for more involved models, bringing into play several 
coupling constants, the resummation procedure is very likely less efficient than for the O(N) models 
since it is performed with respect to one coupling constant only - see below. In this case, the instability 
displayed by the spurious fixed points can be much weaker than for O(N) models. However some of the 
present authors have previously shown [42| that the criterion of (in)stability given above still remains 
reliable in these more general and ambiguous situations. More precisely they have shown, on the 
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example of frustrated magnets (and on the model with cubic anisotropy) that fixed points suspected 
to be spurious from a stability analysis, have been confirmed to be so from additional independent 
arguments. These arguments are: (i) persistence of the fixed point as a non-trivial, non-Gaussian one 
up to (and above) the upper critical dimension d — 4, a fact which is forbidden for a 4 -like theory 
(see [58l [59( 1 and reference therein) , (ii) existence of a topological singularity in the mapping between 
(N, d) and the fixed point coordinates that makes these last quantities multivalued functions of (d, N) 
that is manifestly a pathology. 

From the discussion above it appears that a necessary condition for a fixed point to be considered 
as a genuine fixed point is that it satisfies both the PMS and the PFAC. We however now show on 
the example of the 2d O(N) models that, although necessary, these conditions are not sufficient. 



C. O(N) models in two dimensions: anomalous apparent convergence 



The same kind of analysis of the perturbative results obtained from the (/> 4 model in three dimensions 
can be performed for all N in two dimensions. We show in FigUthe exponent w of the two-dimensional 
0(4) model obtained at three, four and five loops in the zero-momentum massive scheme and the 
anomalous dimension 77 at four and five loops (the three loops results docs not lead to a clear stationary 
behaviour). A conformal Borel resummation method has been used. 




5 10 15 20 25 5 10 15 20 



(a) (b) 



FIG. 4: The exponents uj and 77 of the two-dimensional 0(4) model as functions of the resummation parameter 
b at three-, four- and five-loop order (the result for rj at three loops is not displayed since there is no clear 
stationarity for 77 in this case). The parameter a has been fixed at the value obtained from the large order 
behavior: a — a\ a — 0.1789. For ui one has: (a pt,6opt) = (3.1,9) at three loops, (a op t,&opt) = (3.1, 14) at 
four loops and (a op t,&opt) = (3.1,21.5) at five loops. For r\ one has: (a op t,&o P t) = (4.4,10) at four loops and 
(a pt,6opt) = (4.6, 18) at five loops. The dot on each curve corresponds to stationary values of uj — uj(a,b) 
and 77 = uj(a, b) in both directions. 



For this model, because of Mermin- Wagner's theorem [60], the correlation length is infinite at zero 
temperature only and the critical exponents are exactly known: 77 — and ui = 2 [gj |62[ . We can see 
on Fig f4] that although the values obtained for these exponents seem well converged, they are erroneous 
since using both the PMS and PFAC one finds: 77 ~ 0.12 and uj ~ 1.37. It is important to emphasize 
that N = 4 is not an isolated case in this respect. For all two-dimensional O(N) - with N > 1 
models the critical exponents seem to be converged at five-loop order but towards erroneous values. 
For instance, in the Ising model, and at five loops, Orlov and Sokolov 44| have found 7/ = 0.146, and 
Pogorelov and Suslov [63| 77 = 0.145(14) whereas the exact result is 77 = 0.25. 

We have also studied the a-dependence of the critical exponents. We have found here again that 
a op t — ai = 0.1789, see FigfSJ This means that the a-dependence is not either a good indicator of an 
anomalous apparent convergence in this case. 

An analysis of the underlying reasons of this anomalous convergence has been performed in the 
Ising case in [45[ (see also [gj, Ha|)- The explanation is that there very likely exist, in the 0- function, 
terms such as 1 — (1 — u/u*) e with u* the fixed point value of u and e a small number (probably 1/7 
in the Ising case) |45l |66| . In the perturbative expansion performed around u = 0, such terms lead to 
small contributions to the /3-function that seems to be under control. However they play an important 
role in the vicinity of u* ; they are even non-analytic at this point since their derivatives with respect 
to u are singular at u* . Reconstructing such terms from a perturbative expansion is thus difficult and, 
as a consequence, the perturbative results are doomed to failure although they look converged. Thus 
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FIG. 5: The exponents uj and r\ of the two-dimensional 0(4) model as functions of the resummation parameter 
a, at three, four and five loops (the result for r\ at three loops is not displayed since there is no clear stationarity 
in this case). The vertical line corresponds to a\ ~ 0.1789. The values chosen for a and b are such that the 
exponents are at their stationary value when a = a% (see FigU). Note that the stationary values w(a) and 
77(a) indicated by dots are very close to their values at a = a\ a . 



we are lead to the conclusion that PMS and PFAC are necessary conditions for convergence but are 
not sufficient. 

Let us now make a remark specific to d = 2. In this dimension, the existence of a non trivial 
root u* of the /3 function, stable with respect to the resummation parameters b and a and displaying 
good convergence properties, see FiglHl is not sufficient in itself to know whether the transition is 
trivial (taking place at zero temperature) or not, since u is not directly related to the temperature. In 
principle, the triviality (for TV > 3) or non triviality (for N = I or 2) of the critical exponents should 
be sufficient to conclude. However, as previously emphasized, the presence of strong non-analyticities 
in the two-dimensional j3- functions of the Ising and O(N) models prevent us to do so since they 
completely spoil the determination of the critical exponents. 



u* u 




a 



(a) (b) 

FIG. 6: The fixed point coupling constant u* of the 0(4) model in d = 2: (a) as a function of 6 at fixed 
a = 3.1, (b) as a function of a at fixed b — b opt = 9, 14,22 at three, four and five loops respectively. The 
parameter a has been fixed at a\ = 0.1789. 



III. THE FRUSTRATED O(N) x 0(2) MODELS 

Let us now come to the frustrated models we are directly interested in. A first analysis of the 
convergence of the MS series obtained at five loops in these models was done in [HI]. In the following, 
we study these models by analyzing the perturbative series obtained in the massive zero momentum 
scheme at six loops in three dimensions and at five loops in two dimensions, a case that the MS series 
do not allow to satisfactorily study since the values of the coupling constants at the fixed point are 
out of the region of Borcl-summability. 
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The Hamiltonian relevant for frustrated systems is given by [67H70l j : 

u= f d d a; {i[(a0 1 ) 2 + (a«A 2 ) 2 + m 2 (^ + ^ ) ] + 

(13) 

where (j>i, i — 1,2 are iV-component vector fields. The resummation procedure outlined above can 
be generalized to the case where there are several coupling constants as it is the case for frustrated 
systems. For a function / of the two variables u and v known through its series expansion in powers of 
u and v, the resummation procedure used in (33,111, IzH consists in assuming that / can be considered 
as a function of u and of the ratio z = v/u: 

f(u,z) = Y / an(z)u n (14) 

n 

and in rcsumming with respect to u only. Under this hypothesis the resummed expression associated 
with / reads: 

f R (u, z)=J2 d n (a, a(z),b; z) f°°dt e~* i'^^ (15) 
V Jo [l-oj(ut;z)\ 

with: 



y/l + a(z)u- 1 

u3(u; z) = v — - (16) 
y/l + a(z)u+l 

where, as above, the coefficients d n (a, a(z), b, z) in (|15j) are computed so that the re-expansion of the 
r.h.s. of (|15p in powers of u coincides with that of (fT4"|). Of course, since the resummation is performed 
in only one variable, we cannot expect in this case a convergence of the resummed quantities as good 
as in the O(N) case. 



A. The frustrated models in d — 3 



We recall in Figj7]the results obtained using different approaches. In the (d, N) plane a line N c (d) 
is found in all approaches such that the stable fixed point exists for N > N c (d) and disappears for 
N < N c (d). This result is interpreted as the occurrence of a second order transition for values of N 
above N c (d) and a first order transition for values of N below N c (d). In the e-expansion 35 37[ and 
within the NPRG [H l38W4lT | the lines N c (d) are both monotonic and are very similar (see Figj7]). They 
lead to the fact that N c (d = 3) > 3 and the transition is thus found to be of first order for N = 2 and 
3 in three dimensions. On the contrary, in the MS scheme without e-expansion [28|], the curve N c (d) 
is found to have a S'-shape, see FigfTl and thus, at d — 3, fixed points exist for N = 2 and 3. In the 
massive scheme also fixed points are found for these values of N [3(| . 

The MS scheme perturbative series at five loops were reexamined in p2l |. The bad convergence of 
the resummed series, the analytic properties of the coordinates of the fixed points (u*,v*) considered 
as functions of (d, N) (presence of a topological singularity S, see Fig{7j in the (d, N) plane) and the 
fact that the fixed points found at N = 2 and N = 3 in d = 3 do not become Gaussian when they 
are followed continuously in d up to d = 4 led the authors of [42[ to conclude that these fixed points 
were either spurious or the results non converged. By re-analyzing the resummed series obtained at 
six loops in d = 3 in the massive scheme we show in the following (i) that for sufficiently large values 
of N (typically N > 7) the resummed series for the exponents converge well, (ii) that for N = 2 and 
3 these series do not lead to converged results. The situation is thus similar to what has already been 
obtained in the MS scheme. 



1. The N — 8 frustrated model in d = 3 



Let us start our analysis by the N = 8 case to show how the results obtained at large and small 
values of N for frustrated systems display completely different convergence properties. Since the 
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2.8 3 3.2 3.4 3.6 3.8 4 



FIG. 7: Curves N c (d) obtained within the e-expansion (iV^), the MS scheme without e-expansion (N™) and 
the NPRG approach (iV ( ? ,PRG ). The resummation parameters for the MS curve are a = 0.5, b = 10 and a = 1. 
The part of the curve N™ below S corresponds to a regime of non-Borel-summability. The attractive fixed 
point is called C+ and C+ FD when it is found in the fixed dimension (MS) scheme. 
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FIG. 8: The exponents wi and 0J2 of the three-dimensional frustrated model for N = 8 as functions of the 
resummation parameter 6 at four-, five- and six-loop orders. The dot on each curve corresponds to a stationary 
value of ui = uj(a, b) in both a and b directions with a fixed to its large-order value a\ a — 0.0554. For wi one 
has (a op t,6opt) = (6.5,27.1) at four loops, (a op t,&opt) = (7.9, 16.4) at five loops and (a op t,&o P t) = (7.5, 11.2) 
at six loops. For 0J2 one has (a pt,6o P t) = (5.1, 14.8) at four loops, (a op t,6o P t) = (7.8, 17) at five loops and 
(a pt,6opt) = (7, 12.2) at six loops. 



model involves two coupling constants u and v there are two eigenvalues of the stability matrix of the 
RG flow at the fixed point that we call u± and They represent the generalization of the exponent 
lo of the O(N) models and they rule the stability of the fixed point: it is attractive when uj\ and 102 
have both positive real parts. 

We take for a the value obtained from the large order analysis: a\ a — 0.0554. We find that the PMS 
is satisfied at four-, five- and six-loop orders: for suitable values of the parameters b and a the two 
exponents u\ and o>2 depend weakly on these parameters and are reasonably well converged. This is 
clear from Fig J5] where we show the 6-dependence of ui± and 102- Moreover the difference between the 
values at five and six loops of, for instance, u>2 is small: L02{L = 6) — L02{L = 5) ~ 0.012. Note that 
our values of uj\ and L02 in this case are fully compatible with those obtained in [32[ . 

We have also studied the a-dependence of these quantities and find that the "optimal" value of a is 
close to its large order value (a opt ~ a\ a = 0.0554) as it is the case in the O(N) models, see Fig|9l 

These results indicate that the convergence properties of the N — 8 frustrated model are globally 
similar to those of the O(N) models although less accurate very likely because in the latter case the 
resummation is less efficient due to the presence of two coupling constants. 



11 




FIG. 9: The exponents uj\ and u)2 of the three-dimensional frustrated model for N = 8 as functions of the 
parameter a at four-, five- and six-loop orders. The vertical lines corresponds to the large order value of 
a\ = 0.0554. The values chosen for a and b are such that the exponents are at their stationary value when 
a = a\ (see Fig[8]). 

2. The N = 2 and N — 3 frustrated models in d — 3 

We now analyze the physical values of AT, that is N = 2 and 3. In [42j], it has been found in the 
MS scheme that, because of the presence of the singularity S (see Fig (7]) which exists in this scheme 
for N ~ 7 and d ~ 3.2, the results obtained from the resummed series above and below N ~ 7 are 
very different. In the massive scheme, the series are known in integer dimensions only and it is thus 
not possible to know whether this perturbative scheme leads also to the existence of a singularity. We 
nevertheless show that, within this scheme, the results obtained for N — 2 and 3 are very different 
from those obtained for N > 8 and are fully compatible with those obtained with the MS scheme [42[ . 

Let us first recall that for TV = 2 and N = 3 the fixed point is (in most cases but not systematically, 
in particular for values of a different from 1, 2, 3) an attractive focus, that is u)\ and u>2 are complex 
conjugate and Re(wi)=Re(w2) > 0. We, again, take for a the value obtained from the large order 
analysis: a\ a — 0.1108 for N = 2 and ay Q = 0.095 for N = 3. For these values of a and for L = 4, 5 
and 6, we find that Re(wi) (or equivalently Re(cj2)) considered as a function of a and b is nowhere 
stationary, even approximately, see Fig llOl Moreover, at fixed a and b, the gap between the values 
of Re(wi) at two successive loop-orders: Re(wi)(X + 1)— Re(wi)(L), is always large, of order 0.5 for 
N = 2 and 0.2 in the TV = 3 case, see FigHni Thus neither the PMS nor the PFAC are satisfied for 
these values of N. 

We have also studied the stability of our results with respect to variations of a for characteristic 
values of a and b, see FigfTTl Here again we find no stationarity. 

From these results it clearly appears that the critical exponents deduced from the resummed series 
obtained in the massive scheme in d = 3 display both lack of convergence and of stability for the 
values N — 2 and 3. It seems therefore very likely that the existence of fixed points for these values 
of N is an artefact either of the fixed dimension schemes or of the resummation method. In any 
case, we confirm by this study that there is no reason coming from the fixed dimension approaches 
to question the results obtained either within the e-expansion or the NPRG and that, very probably, 
the transitions found in d — 3 and N — 2 and 3 are always of first order. 

Let us now perform the same kind of analysis for the two-dimensional models. 

B. The frustrated models in d — 2 

As already emphasized, the two-dimensional case is particularly interesting because of the presence 
of topological excitations in the Heisenberg case ■ Because of the homotopy properties of the sym- 
metry group SO(3) of these systems, the topological excitations are different from the 0(2) vortices 
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(a) (b) 

FIG. 10: The (real part of the) exponent uj\ of the three-dimensional frustrated model (a) for N — 2 and (b) 
for N — 3 as a function of b for a = 6 at four, five and six loops. The parameter a has been fixed at the value 
obtained from the large order behavior: a\ a = 0.1108 for N = 2 and aio = 0.095 for N = 3. Other values of 
the parameter a give similar results. 




FIG. 11: The (real part of the) exponent u>i of the three-dimensional frustrated model (a) for N — 2 and (b) 
for N = 3 as a function of a. The vertical lines corresponds to a\ a = 0.1108 for N — 2 and a\ a = 0.095 for 
N = 3. One has taken a — 6 and b — 25, 20 and 15 at six, five and four loops for N = 2 and b = 20, 15 and 
10 at six, five and four loops for N = 3. Other values of the parameters a and & give similar results. 



encountered in the ferromagnetic XY system. It is still an open question to know whether the de- 
confincmcnt of these defects could trigger a genuine phase transition, as in the Kosterlitz-Thouless 
case. Note that such a phenomenon would be surprising since one knows from the spin- wave - low- 
temperature - approach [72H74j that, contrary to the 0(2) model, the spin-spin correlation length of 
0(3) frustrated models is finite at low - but non vanishing - temperature and that vortices tend to 
further disorganize the system. We let aside the delicate question of the very mechanism underlying 
a hypothetical genuine phase transition in these systems and focus on the question of the existence of 
a finite temperature fixed point within the FD formalism. 

As for the XY case the question is to know whether there is a unique or two separate (Ising and KT) 
phase transitions. From the most recent Monte Carlo simulations it has been argued that there are 
two distinct but very close phase transitions, the Ising one taking place at the highest temperature. 
Accordingly one could expect the transition to be characterized by Ising critical exponents. 

In d — 2 and for the values of N > 4 there is no topological defects. As a consequence, there 
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cannot be any other fixed point but the zero temperature one. Thus, for these values of N and 
because of Mermin- Wagner's theorem, the correlation length diverges at zero temperature only and 
with an exponent v which is infinite (exponential divergence of the correlation length). Moreover the 
anomalous dimension 77 is always vanishing at a zero temperature fixed point as can be checked on 
the low temperature expansion performed within the non linear sigma model [72M74j |. Thus, as in the 
two-dimensional O(N) case, any nonvanishing 77 for N > 4 must be considered as an artefact and, for 
the perturbation theory, as a signal of an anomalous apparent convergence, see above. We start our 
analysis by the N = 8 case and then study the physically relevant cases: N = 2, 3. 



1. The N = 8 case in d = 2 

For N = 8 and being given the previous discussion, one should obtain only trivial results as for the 
critical exponents: 77 = and u)\ = u>2 = 2. We have computed 77 and u>i (the largest eigenvalue) as 
functions of the resummation parameters a and b by taking for a the value computed from the large 
order behavior: a\ — 0.0895. We find a stationary solution for the two exponents studied, see Figsfl2l 
However, as in the O(N) models, we find that the value of r\ thus obtained: 77 ~ 0.13 is unphysical 
since it should be zero. We find uj\ ~ 1.79 which is far from the expected physical value u)\ = 2. 





FIG. 12: The exponents u>\ and r\ in the two-dimensional frustrated N — 8 case as a function of b at four- 
and five-loop orders. The dots at each curve corresponds to a stationary value of uj = u(a, b) in both a and 
b directions. For uii one has: (a op t, &opt) = (4.7, 13.3) at four-loops and (a op t, b op t) ~ (4.7, 21.7) at five loops. 
For r) one has: (a op t, &opt) = (4.55, 13.2) at four loops and (a op t, bopt) = (4.55, 22) at five loops. 

We have also studied the a-dependence of these exponents. Here again, we find good convergence 
properties with an extremum around the value a\ a , see Fig ll3l 
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FIG. 13: The exponents u}\ and r] in the two-dimensional frustrated model for N = 8 as functions of a at four 
and five loops. The vertical lines corresponds to a — a opt = 0.0895. The values of a and b are such that the 
exponents are at their stationary values when a — a\ (see Fig|12[). 



It thus appears that there very likely exist in frustrated models, as in O(N) models, nonanalytic 
terms in the /3-functions that spoil the convergence of the resummed perturbative expansion of the 
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critical exponents. We can already assert that this dramatically alters the relevance of the perturbative 
4 -approach for the study of the two-dimensional frustrated systems. 



2. The N = 2 and N — 3 cases ind = 2 



We now perform the same analysis as above for the physically relevant values of N, that is N = 2 
and 3. Let us first notice that, for these values of N, the fixed point starts to exist beyond three 
loop-order only. We fix a at its large order value: a ~ 0.1790 for N — 2 and a ~ 0.1534 for N = 3. 

Let us first discuss the N — 2 case. We find that the correction to scaling exponent uji (and thus 
lu 2 ) is complex for a large range of parameters a and b which means that the fixed point is a focus. 
We show in FigUJ] that there is no value of a and b where lui is stationary with respect to both 
parameters. Moreover, at fixed a and b, the difference between the four- and five-loop results is large. 
This is a clear signal of the nonconvergence of the value of uj\ . In [34J an average value for this critical 
exponent has been proposed: lo\ — 2.05(35) ± iO. 80(55) at five loops. According to our stability and 
convergence principles this value does not really make sense. 
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(a) (b) 

FIG. 14: The (real part of the) exponent cji of the two-dimensional frustrated model for N = 2 as a function of 
b for different values of a (a) at four loops and for a = 1.4, 2.4, 3.4, 4.4 (b) at five loops and for a = 1,2, 3, 4, 5. 
We have chosen a ~ 0.1790. 



The situation is a little bit different for the exponent r\. At four loops r\ is nowhere stationary in the 
a-direction as can be seen on FigfTSJ-a) whereas there is an almost stationary value r\ ~ 0.275 at five 
loops in both a and b directions for a ~ 4.2 and b ~ 11, see Fig[15]-b) (which is compatible with the 
value given in (34j where r\ — 0.28(8)). We have performed the analysis of the stability of our results 
for 77 when a is varied around a\ a at fixed a and b. We find that indeed the five-loop results do not 
vary much with a and that the optimal value of a is close to a\ . 

We conclude that the results for the N = 2 case show no convergence with the loop order and a 
poor stability with respect to variations of a and b but perhaps for the exponent 77 at five loops. Let 
us notice that the value found 77 ~ 0.275 is relatively close to the exact value expected for an Ising 
transition (77 = 0.25). However, at the same time, it is far from the five-loop value rj — 0.146 [44j 
obtained directly with the <t5 4 field theory. We shall develop on this below. 

Let us now examine the N — 3 case, FigUS] The fixed point is again a focus. The difference with 
the N — 2 case is that there now exists a stationary point for Re(wi) at five loops for a ~ 5.95 
and b = 10.25, but not at four loop-order, see Fig ]16I a). where there is no stationarity w.r.t. a. 



At this stationary point one has Re(wi) ~ 1.78 (which is compatible with the result found in [341 ]: 
Re(wi) = 1.55(25) that anyway displays a large error bar). We find stationary points for r\ at four- 
and five-loop orders, see FigfTBT-bV At five loops the value of r\ at the stationary point is 77 = 0.23 
(that compares well with the value 77 = 0.23(5) of [34]). The convergence seems better in this N = 3 
case than in the N = 2 case since now both uj\ and 77 display stationary values. 

From the discussion above one could be tempted to conclude, in the = 3 case, that the value 



i] = 0.23, although affected by a large error bar (drj — 0.05 according to [34J), is sufficiently large to 
ensure that ij does not vanish, as claimed in [34J. In this case the transition would be non-trivial, 
that is, would occur at finite temperature. We now argue that the results obtained at five loops are 
not sufficiently accurate to support this conclusion. The reason is that the error on 77 is, in fact, 
underestimated. To see this we have computed 77(A) (according to our two principles) for all values 
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FIG. 15: The exponent r\ of the two-dimensional frustrated model for N — 2 as a function of b for different 
values of a. (a) at four loops and for a = 0.3, 1.4, 2.4, 3.4, 4.4, 5.2 (b) at five loops and for a = 3.4, 3.8, 4.2, 4.6. 
We have chosen a ~ 0.1790. 




(a) (b) 

FIG. 16: The (real part of the) exponent ui\ and r\ of the two-dimensional frustrated model for N = 3 as 
functions of b at four and five loops, (a) For us\ there is a stationary point only at five loops for (a op t, b op t) = 
(5.95, 10.1) (b) For r\ there are stationary points at four loops at (a opt , 6 op t) = (4.28, 9.6) and at five loops at 
(Qopt, bopt) = (4.45, 16.6). The value of a has been taken equal to its large order value a ~ 0.1534. 



of N between 2 and 8, see FigfTTl As already emphasized there cannot exist nontrivial fixed points 
and, thus, nonvanishing anomalous dimensions rj for any value of N > 4. As seen on Fig |17l this is 
violated by the perturbative results at five loops. This implies that the error Srj on rj at five loops 
is of order rj itself, that is in the N — 4 case, of order 0.20. Being given that rj(N) is monotonically 
decreasing, the error bar increases as N decreases. In the N — 3 case, the error bar is thus at least 
equal to 0.20 and since r](N — 3) is found to be equal to 0.23 it is impossible to conclude that rj is 
nonvanishing in this case. While our considerations extend also very likely to the N = 2 case as for 
the existence of a large error on the result, this case is particular. Indeed for N = 2 one expects 
i] = 0.25 since the transition likely belongs to the Ising universality class. At first sight the value 
found at five loops (77 = 0.275) could seem encouraging. However let us recall that one finds 77 = 0.146 
from the one-component 4 model in d — 2 j44, 63} which is very far from the expected result. As can 
be seen from the > 4 results, the series for frustrated magnets do not exhibit better convergence 
properties than those of the </> 4 model and thus the value of ij found in the N = 2 frustrated case 
should very likely be interpreted as a numerical coincidence. This conclusion is reinforced by the fact 
that, as explained previously, the stability properties of 77 in the N = 2 case are also unsatisfactory, 
see FigHni 



IV. CONCLUSION 



We have investigated the series obtained from FD perturbative approaches of O(N) models and 
frustrated magnets both in d = 2 and in d = 3 at five- and six-loop orders respectively. From a general 
point of view the result of our study is that only the O(N) models in d = 3 provides unambiguous and 
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FIG. 17: Exponent 77 of the two-dimensional frustrated model as a function of the number of spin components 
N. 



precise results. For frustrated magnets, our results in d = 3, that show an absence of stationarity of 
the exponents considered as functions of the resummation parameters a and b and a bad convergence 
with the number of loops, provide strong support to the spurious character of the fixed points found 
for N — 2 and N = 3. Without providing a definitive answer to the question of the nature of the 
phase transition that frustrated magnets undergo in d — 3 our results weaken a lot the predictions 
of a second order behaviour. Since all other studies than the FD approaches (e-expansion, NPRG) 
predict first order behaviour (see [l|) we are naturally led to the conclusion that three-dimensional 
frustrated magnets should always exhibit first order behaviours. 

In d — 2, the situation is more ambiguous since the critical exponents satisfy the PMS in some 
cases (at some orders and for some critical exponents). At first sight, one could deduce from these 
results the existence, for Heisenberg spins, of a finite temperature phase transition triggered by the 
dcconfincment of topological excitations. However a careful comparative study between the O(N) and 
frustrated models shows that the presence of nonanalycities spoils the determination of the critical 
exponents and forbids to conclude. 

It is not clear that only a few higher orders of the perturbative expansion would be sufficient to 
clarify the situation and one has to think about another approach in both d — 2 and d = 3. From this 
point of view the NPRG seems to be able to circumvent the main difficulties. Indeed being not based 
on a perturbative expansion (in the traditional sense) it does not suffer from some of the problems 
encountered in the weak coupling approaches. In particular it seems to be unaffected by the problems 
of nonanalyticities since the value found for rj in the d = 2 Ising case is, within this approach, found 
equal to 0.254 (75|, in excellent agreement with the exact result. The d = 2 case for both frustrated 
Heisenberg and XY systems is under investigation (7(|. 

Let us finally emphasize that the methodology put forward in this article could be relevant for any 
system analyzed within the FD perturbative approach. Indeed, in this case, the existence of spurious 
fixed points is the generic case and one has to be especially careful when fixed points that have no 
counterpart within the e-expansion approach occur. In such circumstances, the principles employed 
here - PMS and PFAC - can be of great interest to reject or to accept the fixed points as physical 
solutions. 
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